Bose-Einstein condensation in inhomogeneous Josephson arrays 
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We show that spatial Bose-Einstein condensation of non-interacting bosons occurs in dimension 
d <2 over discrete structures with inhomogeneous topology and with no need of external confining 
potentials. Josephson junction arrays provide a physical realization of this mechanism. The topo- 
logical origin of the phenomenon may open the way to the engineering of quantum devices based on 
Bose-Einstein condensation. The comb array, which embodies all the relevant features of this effect, 
is studied in detail. 



PACS numbers: 03.75.Fi, 85.25.Cp, 74.80.-g 

The recent impressive experimental demonstration of 
Bose-Einstein Condensation (BEC) |l[ has stimulated a 
new wealth of theoretical work aimed to better under- 
standing its basic mechanisms |^ and, possibly, to exploit 
its consequences for the engineering of quantum devices. 

It is well known that for an ideal gas of Bose par- 
ticles BEC does not occur in dimension d < 2, and an 
"ad hoc " external confining potential is needed to reach 
the required density of states. The same is true for free 
bosons living on regular periodic lattices, while the result 
cannot be extended to more general discrete structures 
lacking translational invariance. 

In the following we shall prove that even for c? < 2 ^] 
non-interacting bosons may lead to Bose-Einstein con- 
densation into a single non-degenerate state, provided 
one resorts to a suitable discrete non-homogeneous sup- 
port structure: indeed, when the bosonic kinetic degrees 
of freedom do not depend on metric features only, the 
particles may feel a sort of effective interaction due to 
topology. The proposed mechanism for BEC in lower di- 
mensional systems is then a pure effect of the structure of 
the ambient space and avoids as well the need of resorting 
to external random potentials as the ones investigated by 
Huang in this is a very desirable feature in view of 
engineering real quantum devices. 

In practice, the behavior of free bosons over generic 
discrete structures is made experimentally accessible 
through the realization of suitable arrays of Josephson 
junctions. The latter are devices that can be engineered 
in such a way as to realize a variety of non-homogeneous 
patterns. We shall show indeed that classical Josephson 
junction arrays arranged in a non-homogeneous geome- 
try - not even necessarily planar - provide an example 
of the proposed mechanism for BEC, leading to a single 
state spatial condensation. 

Theoretical studies of Josephson junction arrays are 
based on the short-range Bose-Hubbard model, since the 
phase diagram of Josephson junction arrays may be de- 
rived 1^] from an Hamiltonian describing bosons with re- 



pulsive interactions over a lattice. In d = 1 the phase 
diagram has been studied by analytical Q and quan- 
tum Monte Carlo methods 0; experimentally, Joseph- 
son junction arrays are used to study interacting bosons 
in one dimension. For a generic array the corresponding 
Hamiltonian is given by 



where A, 



Aij = 1 if the sites 



is the adjacency matrix 
i and j are nearest neighbors and Aij — otherwise; 
a| creates a boson at site i and Ui = a|ai. The phase 
diagram structure reflects the competition between the 
boson kinetic energy (hopping, favouring boson mobil- 
ity) and repulsive interaction (Coulomb, working so as 
to suppress dynamics). In a realistic experimental setup 
Q , the parameters U and V depend on the ratio between 
the intergrain capacitance C and the gate capacitance Co , 
while the parameter J describes Cooper pair hopping. 
Josephson junction arrays allow for a good experimental 
control of C/Cq and J, which can also be varied over 
a wide range. For U ^ J,V , and for bipartite arrays, 
Hamiltonian maps onto the quantum spin-i XXZ 

model 1^. On the other hand, in the weak coupling limit 
("classical " Josephson junctions) \J,V ^ J , realizable 
when C /Co 0, the hopping term dominates the phys- 
ical behaviour of the system, which is then described by 
the tight-binding Hamiltonian 



H = -t^Aija\aj 



(1) 



where t is an effective hopping parameter which accounts 
for a renormalization of the Josephson coupling J. For a 
non translation-invariant geometry of the array the tight- 
binding model, which describes "free " bosons over a reg- 
ular lattice, cannot any longer be interpreted as repre- 
senting non-interacting particles, just due to the ambi- 
ent graph topology. We shall show the dramatic effect 
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of topology already on the simple graph referred to as 
square comb jl^Jl^. This provides an explicit and re- 
markable example of topology-induced mechanism lead- 
ing to a spatial Bose-Einstein condensation in low dimen- 
sion. 

The "square comb " is the graph made of N "fingers " 
of A'^ sites represented in Fig. 1. whose total number 
of sites is 7V^. In the following the generic vertex i is 
labelled with the "coordinate" indices {xi,yi) , i E Zjv 
(where the latter requirement is imposed to guarantee 
periodic boundary conditions). 




Fig. 1. The comb graph: the dots represent Josephson 
junctions and the links describe the topology of the array 
connections, with no reference to their embedding in real 
euclidean space. 



The tight-binding model on the comb graph, is given 
by the Hamiltonian (|^) with the adjacency matrix Aij — 

{Sxi,Xi + l + Sxi,Xi~l) Sy-.0Syjfi + {6y.,y. + l + Sy ■ ^y l) S x^ ,X j ■ 

By exploiting the comb translation invariance in the di- 
rection of the backbone one can perform a Fourier Trans- 
form along the x direction obtaining a new Hamiltonian 
in the variables kx- 



H = -tY^ J2i^y,oSy,y' cos kx + Ayy,)al^yak^ 



■,v' 



(2) 



1 and Ayyi in 



where kx = 2im/N and n = 0,1,-, TV 
the adjacency matrix for a linear chain (i.e. for each 
comb finger). The operator ak,^,y is given by: a^^.y = 
T^^^i^^^"^ '^(x,y)- Notice that (|^) is the sum of N com- 
muting Hamiltonians representing an one-dimensional 
tight-binding model with a local potential at site of 
value: —tcoskx- 

Each of the one-dimensional Hamiltonians appearing 
in eq. (|^) can be diagonalized, in the thermodynamic 
limit. To do this one uses the property that for y > 2 
and y < —2 the eigenvectors of —{SyfiSy^y' cos kx + Ayy>) 



are those of the Hamiltonian describing a free particle 
on the linear chain. Such eigenvectors are: e**'^'"*' with 
eigenvalue — 2cosA:j,, e^*^"*' with eigenvalue — 2coshfcj, 
and (— l)i'e^'^"'' with eigenvalue 2 coshky. If one requires 
that the eigenvalue equations hold also at sites —1, 0,1 
and that the eigenvectors are normalizable, one finds that 
the spectrum of the one dimensional problem is given by 
the isolated point E = — 2isgn(cos kx)^/l + cos^ kx and 
by a continuous part Ro = {E\\E\ < 2t} with a density 
of states given by po{E) = {1/tt){N - l)(4i2 _ £;2)-i/2. 
From the spectra of the one dimensional problems cor- 
responding to different values of kx one obtains in the 
thermodynamic limit the density of states for the tight- 
binding Hamiltonian on the comb-graph (Fig. 2). 




Fig. 2. The density of states of the Hamiltonian (|l|). The 
solid line indicates the continuous part of the spectrum 
po, which is normalized to N{N — 1). The dot lines de- 
note the sets of zero measure, the densities p+ and p_ 
are normalized to N . The x-axis scale is in units of t. 



The spectrum is made of three parts Ra , a — 0, ± : 
for E E Ro the density of states is given by: 

Po{E) = -N{N ~ l)(4t2 - £;2-)-i/2 _ 

TT 

Since J^^^^ pa{E)dE = N{N - 1), in the thermody- 
namic limit, almost all the states, i.e. all the states 
apart a set of measure zero, belong to this region and 
limN^aoi Jj^^^ Po{E)dE)/N^ = 1. In the other two re- 
gions, R^ ^ {E\- VSt < E < -2t} and R+ = {E\2t < 
E < VSt}, the density of states is given by: 



P-{E) ^ p+{E) = p±{E) ^ hN 



\E\ 



V8t2 _ E^VE^ - 4t2 



where h must be chosen so that Jj^ p±{E)dE = 
N. There follows that in the thermodynamic limit 
only a subset of states of measure zero belongs 
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to these regions of the spectrum, in that one has 
limM^oo N-^ Jj^^ p±{E)dE = 0. The states with E £ 
R- play a fundamental role in the study of bosonic par- 
ticles on comb structures. The lowest energy eigenstate 
of (|l|) (corresponding to Eq — —\/8t) is represented in 
Fig. 3. It is constant along the x direction for any fixed 
?/, while along y it decreases exponentially with the dis- 
tance from the backbone. 



A Jro (e/^Cs+tx/s) _ i)V4t2 - E^ ■ 



For T < Tc = iPc)^^ there is Bose-Einstein condensation 
and in Fig. 4 we plot the fraction no of particles in the 
condensate as a function of T for several values of the fill- 
ing (/ = 0.5, 1, 2). The points where the curves intersect 
the T-axis are the critical temperatures for the different 
fillings. 




Fig. 3. The eigenvector corresponding to the lowest 
energy state. It is constant along the backbone direc- 
tion and it decreases exponentially as exp(— arcsh(2) • |a;|) 
along the fingers. 



If one introduces a finite bosonic filling /, i.e. if one 
fills the comb with fN'^ non interacting bosons, fixing 
the number of particles in the gran canonical partition 
function amounts to choosing the fugacity z{N,f3, f) as 



f = N- 



E 



1 



-Re 



- 1 



Po.iE)dE, (3) 



z 



< e 
< cN, 



Since 



with < 

1)"^ p^{E)dE < cN, (c is a number independent of 
N and z), the third term of (||) vanishes in the ther- 
modynamic limit. On the other hand, the first term 
can be positive and finite in the thermodynamic limit 
if z{N,f3,f) g-0tV8 when N oo. Denoting by uq 
the fraction of particles with energy smaller than — 2i, 
no = limjv^oo r^N-^ /^_(0-i(iV)e'3^^ - ir^p^iE)dE, 
one has that hq is a finite fraction of particle condensed in 
a subset of states of measure zero. In the thermodynamic 
limit, if Pc is the inverse temperature for which: 
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dE 



one has that if /? < /3c there always exixts a real and 
positive z solution of the equation 



/ = 



1 



dE 



Ro 



and riQ = 0: there is no condensation. For (3 > jS^ one 
has z = e~^^^ and uq is given by 



Fig. 4. The fraction uq of condensed particles as a 
function of the temperature T for different fillings / = 
0.5,1,2. 



For T — > the order parameter depends linearly 
on T since no oc (T — Tc), and it is analytic in T for 
T close to Tc- The behaviour for T ^ is given by 
no oc VTexp(-(2V2- 2)t/KT). This is different from 
the customary Bose-Einstein condensation in a i-d box, 
where ng = 1 — {T/Tcf"/"^. The critical temperature 
Tc exhibits the dependence on the filling / shown in 
Fig. 5, asymptotically (/ > 1) linear and of the form 

Tc oc — (in/) ^ for / ^ 1. Both behaviours characterize 
a gapped system. The gap (2\/2 — 2)t measures the dif- 
ference between the ground state energy and the bottom 
of the spectral region Rq 
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Fig. 5. The critical temperature Tc as a function of the 
filling /. 



BEC-based encoding and manipulation of quantum infor- 
mation . The idea here is that because of the sponta- 
neous symmetry breaking that characterizes it, a Bose- 
Einstein condensate should be quite naturally described 
by a non- linear quantum mechanics. Such non-linearity 
can be thought of as due just to the effective interac- 
tion of the bosons composing the condensate. In such 
a case the scenario recently described by Abrams and 
Lloyd whereby non-linear quantum mechanics in the 
sense of Weinberg ]l6[ | implies polynomial time efficiency 
in dealing with NP-complete and jJP complex compu- 
tations would hold. It is intriguing that for the system 
presented here the relevant non-linearity, indeed present, 
is due to geometry and topology rather than to physical 
interactions. 



The average energy per particle 



<E>-- 



1 



EdE 



shows that all the particles in the condensate have energy 
— i\/8, namely they are all in the ground state. Summa- 
rizing, in the thermodynamic limit and for T > Tc almost 
all particles have energies between —2t and 2t with the 
distribution {z~'^e^'^ - l)~Vo(£^), while for T < a 
finite fraction no of particles is condensed in the state of 
lowest energy E^ = — VSi- These particles occupy the 
inhomogeneous state described in Fig. 3, in which the 
sites closer to the backbone have larger filling than the 
farther ones. 

In conclusion, we exhibited an explicit example (comb 
array) of Bose-Einstein condensation into a single state 
for non-interacting bosons, induced in dimension d < 2 
by inhomogeneities without disorder and with no confin- 
ing external potential in the discrete geometrical struc- 
ture. The ensuing condensate shows deconfinement in 
one direction, i.e. along the comb backbone, and local- 
ization along the orthogonal direction; this is expected to 
lead to detectable singularities in the response functions. 

The model Hamiltonian used is physically imple- 
mentable by classical Josephson junction arrays, which 
it is possible to engineer in any desired geometric set- 
ting. Thus Bose-Einstein condensates arise as an intrinsic 
device feature, without need of fine-tuning any external 
control parameter. 

Furthermore, it has been evidenced that a comb graph 
may emerge as a relevant geometrical structure if one 
consider a chain of classical Josephson junctions aligned 
along the direction of the backbone and interacting with 
a suitable external environment 

Finally, the same devices, as recently proposed [ p^ , 
might lend themselves to be used for the realization of 
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